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Abstract — In this paper, we address the problem of sta- 
bilization in continuous time linear dynamical systems using 
state feedback when compressive sampling techniques are used 
for state measurement and reconstruction. In [5], we had 
introduced the concept of using li reconstruction technique, 
commonly used in sparse data reconstruction, for state mea- 
surement and estimation in a discrete time linear system. In this 
work, we extend the previous scenario to analyse continuous 
time linear systems. 

We investigate the effect of switching within a set of spar- 
sifiers, introduced in [5], on the stability of a linear plant in 
continuous time settings. Initially, we analyze the problem of 
stabilization in low dimensional systems, following which we 
generalize the results to address the problem of stabilization in 
systems of arbitrary dimensions. 

I. Introduction 

In the last few years, compressive sensing (CS) has 
emerged as a topic of immense interest in the signal process- 
ing, information theory and machine learning communities 
[8], [3]. The overarching principle of CS states that by 
using random, non-adaptive projections, one can accurately 
reconstruct special signals (sparse or compressible) using 
far fewer measurements than what is required using the 
traditional techniques suggested by Nyquist and Shannon. 
Therefore, compressive sensing promises to provide a more 
efficient sensing technique than the existing ones. Reference 
[4] provides an extensive collection of papers related to 
the theory behind compressive sensing, and its applications. 
The single-pixel camera [22], feature-specific imager [1] and 
CMOS separable transform image sensor [21] are examples 
of real sensors that work on principles of compressive 
sampling and sparse acquisition. This paper investigates the 
implications of using such sensors to provide feedback in 
dynamical systems. 

Recently, there has been a growing interest in the control 
community regarding compressive sensing. In [17], the au- 
thors investigate a feed-forward system in which a compres- 
sive sampling system is used to compress the control signals 
using sparse representations. In [18], a networked control 
architecture is improved by using compressive sensing to 
provide more robust packet delivery over a lossy channel. 
In [19], a compressive sensing method was used to improve 
bandwidth limitations on a remote-controlled system. These 
previous papers address utilizing compressive sensing algo- 
rithms to improve pre-existing control architecture. In [11], 
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the authors address the problem of recovering an initial state 
under sparsity constraints. They provide sufficient conditions 
on the number of available observations in order to recover 
the initial state for both deterministic and stochastic systems. 
In a similar vein, the authors in [23] illustrate the technique 
of recovering initial states from sparse observations by using 
a simple diffusion system. 

In this work, we extend the analysis in [5] which in- 
vestigates the use of sensors based on CS techniques for 
providing feedback in control systems. Earlier work [5] 
idealized the notion of sensors using CS techniques, and 
introduced an abstract concept of compressive sensing device 
(CSD). A CSD is a sensor that measures the state of a control 
system. The error in the measurement is dependent on the 
sparsity of the underlying state. Based on the error bounds 
provided by l\ reconstruction techniques [8], the sparsity 
of the underlying state must be within a specified level for 
perfect reconstruction by a CSD. In order to circumvent this 
limitation so that the system performance does not degrade 
in the face of non-sparse states, the paper has proposed 
the use of a sparsifier in order to induce sparsity in the 
states. Furthermore, it has provided a criteria to design a 
linear sparsifier to achieve stabilization in a plant. It was 
shown that for discrete-time linear systems stabilization can 
be achieved if the number of unstable poles is less than the 
level of sparsity that the CSD can handle. The scenario in 
which the condition is violated had remained unaddressed; 
this paper is an effort in that direction. Here, we investigate 
the possibility of stabilizing a plant by switching within a 
class of sparsifiers, each of which fails to stabilize the system 
individually. 

Stabilization of switched systems has been a classical 
control problem which has been studied extensively in the 
past 20 years. Some of the basic problems in stability 
of switched systems have been listed in [16]. Reference 
[15] provides a valuable gateway into the various available 
techniques available for stabilizing switched systems. For 
some linear systems, the approach is to find a Lyapunov 
function which attains a negative value for possible values of 
the states corresponding to each subsystem. Other methods 
include time-dependent switching and using multiple Lya- 
punov functions. Reference [13] considered the case in which 
multiple Lyapunov functions are used, and developed state- 
dependent switching rules for global asymptotic stability. 
Reference [7] provides techniques to stabilize both switched 
and hybrid systems, and [15] contains an exhaustive review 
of the work that has been done regarding the stability of 
switched systems. 



The main contribution of this paper is that we use the 
existing theory in hybrid switched systems to circumvent 
the shortcomings of the system architecture proposed in 
[5], which is based on principles of sparse reconstruction 
techniques. The organization of the paper is as follows. In 
Section 2, we present a brief primer to compressive sensing. 
In Section 3, we present the problem formulation based on 
the architecture proposed in [5]. In Section 4, we provide an 
example of a two dimensional system to illustrate the CSD 
stability concept. In Section 5, we generalize our results on 
stability to systems of arbitrary dimensions. In Section 6, we 
provide simulation results. Finally, we conclude in Section 
7 and provide some future research directions. 

II. Background 

In this section, we present a brief introduction to compres- 
sive sensing techniques and their relation to sparse signal re- 
construction. The content in this section is at times verbatim, 
but a shorter version of Section 2.3 from [14] which provides 
an excellent, comprehensive survey of important results in 
compressive sensing. It is included here as background 
material for the sake of completeness. 

In the CS framework, we acquire a signal xe R N via 
linear measurements 

y = $x + e (1) 

where $ is an M x N measurement matrix modeling the 
sampling system, yG K M is the vector of samples acquired, 
and e is an M x 1 vector that represents measurement errors. 
If x is if-sparse when represented in the sparsity basis , i.e., 
x = ^>a with ||a||o := |supp(a)| < K, then one can acquire 
only M = 0{K\og(N/K)) measurements and still recover 
the signal x [12], [8]. A similar guarantee can be obtained 
for approximately sparse, or compressible, signals. Observe 
that if K is small, then the number of measurements required 
can be significantly smaller than the Shannon-Nyquist rate. 

In [10], Candes and Tao introduced the restricted isometry 
property (RIP) of a matrix $ and established its important 
role in CS. From [10], we have the definition: 

Definition 1. A matrix <E> satisfies the RIP of order K with 
constant S e (0, 1) if 

(l-5)||x|||<||*x||»<(l + *)K (2) 

holds for all x such that ||x||o < K. 

In words, <£> acts as an approximate isometry on the set of 
vectors that are if -sparse in the basis "P. An important result 
is that for any unitary matrix <3>, if we draw a random matrix 
$ whose entries <fiij are independent realizations from a sub- 
Gaussian distribution, then will satisfy the RIP of order 
K with high probability provided that M = 0(Klog(N/K)) 
[2]. In this paper, without any loss of generality, we fix ^ = 
I, the identity matrix, implying that x = a. 

The RIP is a necessary condition if we wish to be able 
to recover all sparse signals x from the measurements y. 
Specifically, if ||a:||o = K, then $ must satisfy the lower 
bound of the RIP of order 2K with 5 < 1 in order to ensure 



that any algorithm can recover x from the measurements y. 
Furthermore, the RIP also suffices to ensure that a variety of 
practical algorithms can successfully recover any sparse or 
compressible signal from noisy measurements. In particular, 
for bounded errors of the form 1 1 e| |a < e, the convex program 

x — argmin ||x||i s.t. ||$x — y ||a < e (3) 
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can recover a sparse or compressible signal x. The following 
theorem, a slight modification of Theorem 1.2 from [9], 
makes this precise by bounding the recovery error of x with 
respect to the measurement noise norm, denoted by e, and 
with respect to the best approximation of x by its largest K 
terms, denoted by x^. 

Suppose that satisfies the RIP of order 2K with 5 < 
V2 — 1. Given measurements of the form y = Q^x + e, 
where ||e|| < e, the solution to (3) obeys 

||x-x|| 2 <C e + Ci l|X ~-* K|1 (4) 

where 

n _ 4(1 + 3) n _ l + (V2-l)a 

l-(V2+l)d' 1-(V2 + 1)5 

While convex optimization techniques like (3) constitute 
a powerful method for CS signal recovery, there also exist a 
variety of alternative algorithms such as CoSaMP [20] and 
iterative hard thresholding (IHT) [6] that are known to satisfy 
similar guarantees under slightly stronger assumptions on the 
RIP constants. In this work, we assume that measurement 
noise is absent, i.e., e = 0. 

III. Problem Formulation 

In this section, initially we introduce the control archi- 
tecture as proposed in [5]. Based on this architecture, we 
formulate the problem statement. 

As introduced in [5], the CSD is an idealized device 
that works on the principles of li reconstruction algorithm 
[8]. It is assumed that the reconstruction error of the CSD 
obeys Q. Therefore, if the input state is S'-sparse, there is 
no measurement error i.e, the output of the CSD perfectly 
replicates the input. Figure [T] illustrates the architecture of 
the control loop when a CSD is used to provide linear state 
feedback to stabilize a plant. In [5], the idea of using a 
sparsifier, G, has been proposed, which replaces some of the 
state measurements with zeroes instead of their true value. 
This increases the sparsity of the resultant state vector to the 
desired value so as to ensure perfect measurement using the 
CSD. However, we can observe that if the number of unstable 
poles are greater than the maximum allowable sparsity, there 
are instances for which stabilization cannot be ensured by 
linear state feedback law. An example of such an instance is 
as follows. Consider the following decoupled 2-dimensional 
system: 
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Fig. 1: CSD Dynamical System Framework 
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Let us choose the sparsifier G = 

[ki k 2 ; k 3 fc 4 ] be the feedback gain. The closed-loop 
system then has the following form: 
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from which it is obvious that the pole corresponding to A = 2 
cannot be shifted for any value of K. Therefore, the system 
cannot be stabilized by providing linear state feedback by us- 
ing the above sparsifier. This leads to the following question: 
Is it possible to construct a class of sparsifiers, independent 
of the structure of the system matrices, A and B, for which it 
might be possible to stabilize the system by switching within 
the class? We assume that any single sparsifier in the class 
leads to an unstable closed-loop system. In the following, we 
construct a class of sparsifiers that depends on the maximum 
sparsity in the state that can be tolerated by the CSD for 
perfect reconstruction, denoted by S, and the dimension of 
the system. 

Before defining the class of sparsifiers Qs, let us introduce 
a vector shift operation which we denote by rshift(a,b), 
where a E R™ and b € Z+. rshift(a,b) performs a right 
shift operation on each element in the vector a and substitutes 
zeroes in the shifted slots. Let a — [ai, da, ci3, • • • , a n ], then 
rshift(a,b) = [0, . . . , 0, a 1+b , a 2 +b, ■ ■ ■ , a n -b] where the 
output of the operation is still a vector in K™. 

We now define a class of sparsifiers that introduces n — S 
zeroes in the state measurement, where n is the dimension 
of the state space. Let the set of sparsifiers be denoted by 
Qs = {G l }, where i G {1,2,..., [|]} and d e R nxn . For 
each i £ {1,2, [%~\ — 1}, define each Gi as a diagonal 
matrix with the diagonal denoted as T>i £ K™. T>i is defined 
as follows 



9,£K ifj={l,2,...,L|J} 
else 



Forze {2,3,..., Tgl }, V i = rshift(V 1 ,ix [fj). 

As an example, let us consider the case when n — 5 and 
S = 2. The set Qs then consists of the following set of 
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Each sparsifier in the set Q g introduces at least 3 (n — S) 
zeroes in the state. 

Now we present the problem statement. Consider the 
control architecture shown in Figure [2] The pair (A, B) is 
assumed to be stabilizable. Moreover, we assume that there 
is no K such that A c i = A — BKGi is stable for any 
single Gi € Qs- We want to find a switching scheme among 
the sparsifiers in the class Qs such that the system can be 
stabilized. 

Fig. 2: Switched System with CSD 
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In the next section, we investigate the problem for a two 
dimensional system. 

IV. Two-dimensional System Example 

Let us consider a two dimensional system with sparsified 
feedback. In this section, we assume B = I for the sake 
of simplicity. We define two sparsifiers in this example, and 
show a stability result which motivates our main theorems 
in the next section. 

Let 5s be defined as follows: 
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The closed-loop behaviour of the system is given by the 
following equation. 

x = Ax + Iu = Ax — Kx 

= Ax- KG l x = (A — KGi)x, i e (1, 2) 

From [15], we know that there exists a switching scheme 
to guarantee stability if there exists a convex Hurwitz combi- 
nation of the subsystems. We denote this convex combination 



by A. 



A = ctiAi where ai = 1 



(6) 



Next, we investigate the existence of a convex combination 
for a general state matrix A of the following form. 



A = 

Let K = kl, where k £ M. 

A a = 



a\ — k a,2 
a 3 a 4 



ai a 2 
a 3 a 4 



A 2 = 



ai a 2 

03 Cl4 — fc 



Let us choose a, = ~ Vi. Therefore, from ^ we obtain 
the following value of A. 
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Let A be an eigenvalue of A, and v be the corresponding 
eigenvector. 



Av = Xv 



K—lV = —V 

2 2 

n 2 
3t; = (A-~)w 



(7) 



From the above equation, we can conclude that the eigen- 
values of A can be arbitrarily placed by choosing a proper 
k. For k > 2Re(X max (A)), we obtain a Hurwitz convex 
combination of the two subsystems. 

In the next section, we generalize the result to systems of 
arbitrary dimensions. 



V. Proof of Stability 

Consider the system x = Aj,x,i £ (l,2,...,n). Also let 
Ai be non-Hurwitz Vi. 

Theorem 1: If the matrices Ai have a Hurwitz convex 
combination, then there exists a state-dependent switching 
strategy that makes the switched linear system quadratically 
stable [15]. 

Proof: If there 3 a convex combination such that 
A = Y^i=i a iAi where 2~27=i &i = 1 is Hurwitz, then by 
the definition of Lyapunov stability, for each positive definite 
matrix Q, there exists a positive definite matrix P such that 



-Q 



This implies that 

n n 

C^aiAifP + PC^aiAi) 

i=l i=l 
n 

J2 a i( A T P + PA i) = -Q 
i=l 

n 

=> x T ^2ai(AjP + PAi)x = -x T Qx < 

Vx e M n \ {0} 

=> 3i £ [l,n] s.t. x T ai(AjP + PAi)x < 
Vx £ R n \ {0} 

Let Q{ be defined as follows: 

a i :={x:x T {AjP + PA i )x<0}, i=[l,n] 

Let x £ 7Z" \ {0}. Then for a given positive definite matrix 
Q, -x T Qx < 0. Since there exists a convex combination of 
A = 57J"_j ctiAi which is Hurwitz and satisfies the following 
Lyapunov equation 

-Q = A T P + PA 

n 

=> -Q = A T P + PA = ^a l {AjP + PA l ) 

i=l 

n 
i=l 

=4> 3 at least one i G 7? s.t. x T (AfP + PA t )x < 

=*> a; G U ieV fli 

On the other hand, it can trivially be shown that x £ 
Uj e [i )n ]Qj x £ lZ n \ {0}. Therefore, the following 

holds true: 
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The switching strategy is based on keeping the system x — 
AiX active in fij since it will decrease the function V(x) := 
x T Px along solutions to the differential equation. ■ 
Our next goal is to show that for linear stabilizable dy- 
namical systems, the class of sparsifiers Qs can always make 
the convex combination of subsystems in a CSD dynamical 
system framework Hurwitz. Before we start, let us define the 
following: 



Ai = A-KG h i £ [!,.,[- 



(8) 
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Theorem 2: Given a stabilizable dynamical system and 
the class of sparsifiers Qs, there always exists a convex 
combination of Aj's, which is Hurwitz. 



Proof: The convex combination is defined as A — 
J2i otiAi where £\ = 1 

3 - j2 a i A i=Yl a ^ A - B&G ^ (9) 

= ^OiiA-BKGi) 



A-BKY, 
diag(ai5i, . 



,ai 5 i...a rf1 <7 rs1/ 



Since the system is stabilizable, 3K such that(A — BK) is 
stable. Moreover, if we choose <Xj > and > Vi then 
IT 1 exists. Therefore, K = KT,' 1 . ■ 
Based on Theorem 1, it is clear that there exists a 
stabilizing switching control scheme because there exists a 
Hurwitz convex combination of the individual subsystems. 

VI. Simulation Results 

We consider the following system matrices, and the spar- 
sifter class. 
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Clearly the system (A, B) is stabilizable (in fact, control- 
lable). A linear controller based on the gain 



K 



will make A — BK Hurwitz. Without loss of generality, we 
let K = I, Letting <Xj = |, Vi, and g\ = g 2 = 8 leads us to 
the following: 

U o" 



KY, = 



4 



By Theorem 2, we can find a class of sparsifiers Qg that 
yields a Hurwitz convex combination of the subsystems. 
Using Theorem 1, we can find a common Lyapunov function, 
and a state-dependent switching strategy. 

" " -3 



Let Q = I. Substituting A = 



-2 



in the Lyapunov 



equation A P 
P: 



PA = — Q leads to the following value of 
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This creates two conic regions which determine the switching 
strategy between the sparsifying transforms. 
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Figure [3]is the state trajectory with initial value of xq = [2 
1] T . Figure |4] shows the phase portrait. The portrait shows 



Fig. 3: State Trajectory 

State Evolution 




Fig. 4: Phase Portrait 
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the trajectories moving to the origin in a region where both 
subsystems decrease the Lyapunov function. 

Figures [5] and [6] show the state evolution and phase 
portraints, respectively, for B — [1 1]'. The initial 
condition is chosen as xq = [2 -l] . Once again, 
stability is observed via the state-dependent switching 
strategy. 



Next, we consider the following coupled linear system: 
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The system has both eigenvalues at A = 1. The system is 
stabilizable (in fact, controllable). A linear controller based 
on the gain 

K = [8 8] 



Fig. 5: State Evolution 
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Fig. 6: Phase Portrait 
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will make A — BK Hurwitz. Repeating, as before, we 
conclude that if ajj = | Vi, gi = g 2 = 8, and K = [l l] 
we obtain the following equation: 

= [8 8] 

Following the same procedure as before, we obtain the 
following value for P. 



P 



6.5 
1.75 



1.75 
0.75 



n 2 



(x : x T {AlP + PA x )x < 0) 
(x : 8x1X2 + 5.T2 — 15.Xj < 0) 
(x : x T {AjP + PA 2 )x < 0) 
(x : 13x1 — 8x1X2 — 1x\ < 0) 



P generates two conic regions which determine the switching 
strategy between the sparsifiers. Figures [7] and [8] show the 



state evolution and phase portraits respectively. The initial 
values is chosen to be xq = [1 — 0.5] T . 

Fig. 7: State Evolution 
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VII. Conclusion 

In this paper, we have investigated the effect of switching 
within a class of sparsifiers, introduced in [5], on the stability 
of a linear plant in continuous time settings. Initially, we 
analyzed the problem of stabilization in low dimensional 
systems; based on which we addressed the problem of 
stabilization in systems of arbitrary dimensions. A key 
contribution in this paper is the construction of a general 
class of linear non-invertible sparsifiers, and a corresponding 
switching strategy to stabilize the overall system. 

Some of the future directions for research include finding 
other classes of sparsifiers which might be useful in stabi- 
lizing both linear and non-linear dynamical systems. We are 



currently addressing the problem of limiting the number of 
switches. We believe that a cost on the number of switches 
will provide a metric on different classes of sparsifiers, and 
shed some light on the problem of finding the optimal class. 
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